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LIQUID BEHAVIOR IN THE RESERVOIR OF THE 
SOUND SUPPRESSOR SYSTEM 

SUMMARY 


A large water reservoir of rectangular plan form is part of a sound 
suppression system projected for use in the static firings of the booster stage 
of the Saturn V launch vehicle. Upon ignition of the rocket engines, a surge 
wave is created in the water mass at the discharge side of this reservoir. The 
surge wave will travel across the water surface and be reflected on the opposite 
wall, and thus travel back and forth until it is eventually damped out by wall 
and internal friction. 

The characteristics of the surge wave are investigated by treating 
first the response of the liquid to harmonic excitation of one' side wall of the 
reservoir; "admittance" is thus determined, with which the response of the 
liquid to any arbitrary motion of the container wall can be obtained. Motion is 
investigated for a rectangular and a sinusoidal velocity respectively of one 
container wall. The liquid is considered incompressible, nonviscous, and 
irrotational. 

The velocity potential, free liquid surface displacement, velocity and 
pressure distribution, fluid force and moment have been obtained for various 
pulse durations; these results serve as design values for the reservoir struc- 
ture. The wave form was computed electronically for a time coordinate with 
one-sixteenth second increments. These increments yield the wave motion in 
true time by assembling them in a film which shows quite lucidly the behavior 
of the liquid surface. 


SECTION I. INTRODUCTION 

A reservoir of five-million gallon capacity is part of the projected 
sound suppression system of the Saturn V Static Test Facility. By injecting 
large quantities of water into the exhaust jets near the engine nozzle exit, the 
sound power generated by the exhaust jets will be reduced considerably (Fig. i) . 
The high velocity of the exhaust creates a suction which causes the water to 
flow, but also creates a surge wave in the container when the rocket engines 



are ignited. The water system reacts as if one wall were suddenly moved. 
However, the water level is kept constant during the operation; the surge wave 
is of two-dimensional form. For this reason the motion of the liquid in the 
rectangular container can be treated as that in a rectangular container of in- 
finite width. 

The container is considered to be of length a and filled with incom- 
pressible liquid to a height h. To describe the real physical effect of the suc- 
tion created by the firing of the rocket engines, one wall of the container is 
assumed suddenly to move for a finite duration. The displacement y 0 of the 
wall caused by this motion is considered small compared to the length of the 
container; the level of the undisturbed liquid thus remains in the average at its 
location x = 0 (Fig. 2). 

The problem posed is to determine the characteristics of the surge 
wave. To resolve this problem, the characteristics of the liquid motion are 
investigated by treating first the response of the liquid to harmonic excitation 
of one side wall of the basin. By determining the "admittance, " the response 
of the fluid to an arbitrary motion of the container wall can be obtained. This 
is performed for a rectangular and a sinusoidal velocity respectively of one 
container side wall. The velocity potential of the liquid is determined both 
for the time interval of the pulse and for the time after completion of the pulse. 
The free fluid surface displacement, velocity and pressure distribution, fluid 
force and moment can then be determined and serve as design values for the 
reservoir structure. 

The motion of an ideal incompressible liquid in a rectangular container 
of infinite length filled to a height h is then considered. A cartesian coordinate 
system Oxyz is introduced such that the yz plane is located in the quiescent sur- 
face of the liquid perpendicular to the gravity vector. The x-axis is pointing 
vertically upward. The height of the free liquid surface above the plane x = 0 
is called x = x (y,t) ; and is caused by some disturbance of the container or 
one of its parts, such as a side wall. 

The velocity of the ideal, incompressible liquid in the container can be 
represented as a gradient of a velocity potential tf> , which, due to the continuity 
equation has to satisfy the Laplace equation. The normal velocity of the liquid 
at the container walls is equal to that of the container wall. The free fluid sur- 
face boundary condition is described by the kinematic and dynamic condition, 
of which the first expresses that the normal velocity of the fluid particles at 
the free surface is equal to the normal velocity of the free fluid surface. The 
latter condition is obtained from the unsteady Bernoulli equation for the pres- 
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sure p equal to the ambient gas pressure p Q . It is the first integral of the 
Euler equation. 

The problem is then to find those solutions from the class of harmonic 
functions satisfying the container wall boundary conditions which also satisfy 
the kinematic and dynamic condition. By proper transformation, the wall 
boundary conditions can be made homogeneous, and the problem is then to 
solve the Poisson equation with these boundary conditions. 

The problem of free fluid oscillations in a rectangular container has 
been treated by Rayleigh [1] , in 1876.. In recent years, the problem of forced 
fluid oscillations has grown in importance [2, 3, 4] . Lorell [5] gave the 
linearized flow of a liquid in a rectangular container of infinite length for trans- 
latory harmonic excitation of the complete container. The purpose of the 
present investigation is to determine the liquid response (such as free fluid 
surface amplitude, pressure and velocity distribution, fluid force and moment) 
caused by pulse excitation of one container side wall. 

SECTION II. RESPONSE OF LIQUID IN A RECTANGULAR 
CONTAINER OF INFINITE WIDTH DUE TO HARMONIC 
EXCITATION OF ONE SIDEWALL 

Before proceeding to. arbitrary wall excitation the solution for harmonic 
excitation must be known. Treating the liquid as incompressible, irrotational 
and inviscid results in a presentation of the velocity as a gradient of a velocity 
potential <fi that satisfies the Laplace equation 

V 2 <p = 0 . (2.1) 

This equation has to be solved together with the boundary conditions of the pro- 
blem. At the tank walls the normal velocity of the liquid and the container are 
equal to each other; after linearization it is obtained 


o 

II 

21 £ 

at the side wall y = 0 

(2.2) 

9 <b ._ ifit 

= y ifie 
9y o 

at the side wall y = a 

(2.3) 

|4 = o 

9x 

at the container bottom x = -h . 

(2.4) 
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The free fluid surface condition is obtained from the linearized kine- 
matic condition 90/ 9x = 9x/ 9t and from the dynamic condition for p = 0 ob- 
tained by eliminating the free fluid surface displacement x from the linearized 
instationary Bernoulli equation 90/ 9t + gx + p/p = 0. Thus, 

+ S = 0 at the free fluid surface x = 0. (2. 5) 

Here, g is the gravity constant. The displacement y^ is considered small in 

order to maintain the undisturbed free fluid surface in the average at x = 0 
(i.e. , a.h = (a ±y ) • h* which yields h' 1 ' « h(l+y 0 ) ; therefore, it should 
° IT 

be y / a « i) . 


The transformation 
0 (y, x, t) = |*(y,x) + ifiy^ 


v^ i o t 

r,x,t)= *(y,x) + e 


( 2 . 6 ) 


makes the boundary conditions (2.2) and (2.3) homogeneous. Instead of the 
Laplace equation (2.1), one has to solve the Poisson equation 
ifiy 

V 2 * = - — — 5 (2. 7) 


a 


with the transformed boundary conditions 


9i£ 

9y 

9 £ 


9x 


g 


9x 


0 

0 

- = ma y 0 i 


at y = 0, a 
at x = -h 


at x = 0 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 


A solution satisfying the boundary conditions (2.8) is given by 



Introducing this expression into Poisson's Equation (2.7) yields 


A "(x) - — A (x) = 0 for n = 0 (2.12) 

n . ar n y ' 

and 

y 

A "(x) = -iS2— . (2.13) 

o a ' ’ 

This infinite set of ordinary differential equations in x has to be solved 
with the boundary conditions (2.9) and (2. 10) which with (2.11) yield 


A ' ( -h) = 0 for n = 0, 1, 2, (2.14) 

iS2 3 ay 

gA n '(°) - £2 2 A n (o) = 2(-l) n - 2jr 2 ° for n = 1, 2, (2.15) 

iS2 3 y a 

gA Q '(o) - S2 2 A q (o)= g 2 - . (2.16) 


Here, the right-hand side of equation (2. 10) has been expanded into 
a Cosine-Fourier-series and terms of the left-hand side of equation (2. 10) have 
been set equal to those of the right-hand side. These become 


i£2 


,.2 

3 y -2- = 

y o 2a 




where 


a 

o 


i£2 3 ay 

l_o 

3 


2a(-l) n iS2 3 y o 

and oi = o — o 

n 7i rr 


(2.17) 


( 2 . 18 ) 


The solution of equation (2. 12) with the boundary conditions (2. 14) 
and (2.15) together with r] = Si / io as the ratio of the exciting frequency to 

the natural frequency of the liquid and co ^ tanh ( n 7 r — j ( see Table 1 . ) 


is 


A n (x) =2(-l) n 


iftay^rj 2 ^ cosh (x + h) 

7r 2 n 2 ( 1 - r \ n 2 ) cosh j 


(2.19) 
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while that of equation (2. 13) with the boundary conditions (2. 14) and (2. 16) 
is 


y iy , ift ay 

A o <x > = - in ^ (x2+ 2xh) ■ ir ■ (2 - 20) 

i£2 1 

The velocity potential <p (y,x, t) caused by the harmonic excitation y^e ol 
one tank wall (y = a) is, therefore, with (2.6) given by 


i £1 1 

0(y,x,t)=y Q e ift 



a _ gh 
6 aQ ‘ 


2a 


( x 2 + 2hx ) 


2a v / v n n 

1 <-» ITT 

n=l 


cosh [ J a L (x + h)] 


V ) 2 u/ h \ 

n n cosh(n7r — ) 
a 


cos 



( 2 . 21 ) 


The free surface displacement, pressure distribution and velocity 
distribution, as well as fluid force and moment can be determined from the 
potential by differentiations and integrations with respect to the time and spatial 
coordinates. 


The surface displacement of the liquid which is measured from the 
undisturbed position of the fluid is (Figs. 3, 4, 5 and 6) 


x (y, t) = 



00 



( - 1)% n 2 



gh_ 

aft 2 


( 2 . 22 ) 


For very small ft « 1 and very large ft » 1, the liquid surface regresses like 

- — v t The pressure distribution at a depth ( -x) is given by 
a o 


P = 


3 ch 

- p -rf -pgx = - pgx + py 


£ 

2a 


2 iflt r 

ft 2 e 

° L 

oo n 2 

2a y . n _ n § 

+ x2 „=l n* cosh (jaA) 


a gh 1 


ft 2 a 


2a 


(x 2 + 2hx) 


cosh [ja 21 ( x + h) J (m 
cos — y 

2 I — !L 1 V a 


')] 


(2.23) 
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The pressure distribution at the wall y = a is, therefore, 

1 


P _ = - Pgx + Py n 2 e 
y— d. u 


_a_ _ gh_ 

3 fi 2 a 


2 a 


(x* + 2hx ) 


2a 
7 r 2 


V cc 

(1 -77 2 ) ; 

n=l ' n n 


cosh La 31 (x + h)] 


2 cosh 


(f) 


and at the wall y = 0, it is (Figs. 7, 8 and 9) 

a gh J_ 


_ 2 lflt 

P = - pgx + p y R e 
y=0 o^o 


2a y> 'n 

* 2 L (i -o 


00 ( - 1 ) n T7 , 2 


6 £2 ‘‘a 2a 

cosh [^r ( x + h) 

? , f mr h ^ 

n cosh I — - — J 


( x 2 + 2hx ) 


At the container bottom (x = -h) the pressure distribution is 


(2.24) 


(2.25) 


bottom 


pgh + py Q 

« mt 

£2 2 e 

oo 

(-D \ 2 

2a v 

+ _ 2 L 


n=l 

(i-r7 n 2 


yl a gh hi 


2a 


2 \ „2 


STa 2a 


(^)J 


(2.26) 


The fluid force is obtained from integration of the pressure along the sidewalls 
o 


F = / (p - P \ 

y \ y=* y=°) 


dx 


(2. 27) 


which with the fluid mass m = p ah yields the expression (Fig. 10) 


2 1 

F = my Q,^e 
v o 


1 4 tanh [ ( 2n - 1 ) tt a ) 1 ^2n - 1 


— + 


2 


n=l (2n - l) 3 tt | ^1 - v 2 2n _ t 


(2. 28) 


The moment of the liquid about the center of gravity of the undisturbed 
fluid is given by 
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M = 
z 


= / (P y=a - P y=0 ) (| + x ) 0* - h bottom ( y - I) 


and yields the expression 


M = mafl 2 v e 


iS2t 


(2n - 1) 

The velocity distribution is given by 



OO 

fa 

2 

1 


24h/a 

tH 

1 

1 

Zls 

' 


1 

2 

A 

-D^j 

75 : = 

2 Vcosh 

^2n - 1) n il ^ 

- 1 

J 


— (2n - l) 4 


d(h it2t 

v = 7T - y if2e 
9y -'o 


(2.29) 


I _ 1 V , n n coshfr (x + hxK fa 

a 7r (w h \ ( 2 \ la 

n oosh (“) M (2.3ft) 


n=i 


At the left container wall (y = o) the velocity v is zero and at the right con- 
tainer wall (y = a) it is y^ifie 1 ^* , which agree with the given boundary con- 
ditions . The vertical velocity distribution is 


u 


_ 9<i> _ 

ax 


y ifie 
o 


i£2t 


x+h 


+ 7 2 <-» 

n=i 


n sinh (x+h)] 77 !! 


n cosh 


y ' ' x / cos 

P)(‘ - 0 


fr 1 ) 


(2.31) 


At the container bottom (x = -h) the velocity is u = 0, which agrees with the 
boundary condition. 

SECTION III. RESPONSE OF LIQUID DUE TO 
ARBITRARY EXCITATION 

To determine the response of a liquid with a free fluid surface in a 
rectangular container of infinite length with one moving side wall caused by 
an arbitrary nonperiodic excitation y(t) , the function y is written into its 
simple harmonic components by means of the Fourier integral 
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The solution to the steady state problem yields the "admittance" of the fluid 
system, which is combined with the components of the response of the liquid 
to obtain the response of the system to the original nonperiodic excitation 
functions . Then 


1 °° °° 

$(y,x,t) = ~ f f $ (y> x; £2 ) y(r) cos£2(T-t) 

TT J •’ 


dTd£2 


o -00 


is the velocity potential of the liquid due to the nonperiodic excitation y(t) . 
Therefore, by rewriting the potential (2.21) as 


cj) ( y, x, t) = y^cos Qt 


+ 


2a 


yL 

2a 

a 

6 

(x 2 +2hx) - 

gh 

£2 2 a 

(-1) 

n 2 
^n 

cosh [tF" (x + h)] 

fW ^ 

(1 - 

2\ 

^n 

n 2 cosh ( ] 

COS yj 



V a ) 

_ 


the n admittance n <p is given by 


y 2 _ a 1 

2a 6 2a 


<$> (y,x;£ 2 ) = ir - 7 - TZ (x 2 + 2hx ) - 


gh 

£2 2 a 


2a 

7T 


v , .,>n £2 2 cosh [a 21 (x +h)J fm \ 

* c„ sh (as.) oos M 


and the velocity potential for an arbitrary nonperiodic excitation y(t) yields 

00 +00 


<My,x,t) = ^ ff 


0 —00 


Y — — _A_ / ^-2 I or Y \ 

2a 6 2a U 2hx) £2 ’a 


2a n Q 2 cosh f¥( x+h)J fmr v\l . 

cos MJ s 


y(r) cos £2 (r-t) dTd£2 , 



A, RESPONSE OF LIQUID DUE TO SINGLE RECTANGULAR PULSE 


A single rectangular pulse of magnitude and duration t 
ing at t = o (Fig. 11) can be described by the Fourier-Integral 


OO CO 

y(t) = v(t)=- f [ v(r ) cos 0 (t - r) dr 

IT J J 


dfi 


O -oo 


With 


S’ 


V (t) = < 


for 


for 


o < t < t 


t > t 


V. 

the expression obtained is 

V °o 


J 


, L . o r cos fitsinfit* , or 

(t) = — I d fi + — I 

7T J £2 7T J 


V oo 

£ r sin fit ( 1 - cosfitp) 


fi 


dfi 


The velocity potential for rectangular pulse excitation is, therefore, 

= Zo _ a _ ± 


$ 


(y,x,t) = -7 / 

7T J 


|_ 2 a ” 6 2 a 


. o .1 cosfitsinfit 0 
( x 2 + 2 hx) dfi 

J ui 5 


v <=° , sinfi t cos fit 


or h 
7 r J ** a 


2 av oo 


fi- 


dfi 


7T 


v °°r 2 a "i 

- f -? - ! - T- (x 2 + 2 hx) 
ir J [ 2 a 6 2 a 


o v . n coshta?" (x + h)]cos (tt y) r 

! M' 11 „ J 


oo fi cos fit sin fit 


n 2 cosh { — ^ j 


(w " - fi") 
o n 


sin fit £1 - cosfitp J 


fi 


dfi 


o° sin fit £ 1 - cosfitp J 


dfi 


v r gn <- 

tt ^ a fi ' 

o v> , . , n cosnL ( x + h)J cos (tF y) r fi sin fi t Cl — cos fi 

+ -TT- L (- 1 ) : >n/h -i J 


o 

2 av_ 

TT 


n cosh[ 


n 2 cosh 


~ rm_ h ~ 
_ a _ 


(oj " - fi") 
n 


start- 


dfi 


( 3 . 1 ) 
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methods 

oo 

/ 

O 

oo 

I 

0 
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O 
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O 
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I 

O 
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The occurring integrals can be solved with complex integration 
( see Appendix) . 


The solutions are: 


cos £2t sin £2t 


£2 


sin £2t cos £2 t 
c 
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d£2 


d£2 


sin £2t 
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sin £2 1 £ 1 - cos £2 1 J 
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sin £2 1 cos £2 1 
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£2 cos £2t sin £2t 
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r 


v. 
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7T 

2 
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t < t 

o 

o 

for 

t > t 

o 

o 
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t <t 

o 

£ 

2 

for 

t > t 

o 

7T 



2 

for 

t > o 

£ 

2 
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t < o 

7 t ( 2t -t) 
4 o 
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I t2 

4 o 
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4<‘ 2 

+ t 2 ) 
o 

for 

f- 

t 

o 

for 
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- — cos co t cos co t for t < t 
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n o 


7T 

— sin co t sin co t for t > t 
2 n n o c 
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°° fi sin fit cos fit 
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co 2 - fi 2 
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dfi 


fi sm fit 
— 2 — 7ZT~ dfi 
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t sin 
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o 
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t cos 
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oj 
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cos CO 

t 

for t > o 
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With these results, the velocity potential for excitation caused by a 
single rectangular pulse of the container wall y = a is given by: 


4>(y,x,t) = v 




yl _ a _ 

2a 6 2a 


_ , , 2t 2 - 2tt + t 2 

(x z + 2hx)l + ^ ^ ° ° 


* - 2 


J * \ 2tt -t 2 
o o 


(3.2) 


n=l 


. rmr . , , .1 (rm \ (cosco (t -t) + cosco t] 

n cosh \_~ar (x + h)J cos \1r y) I no n I 


n 2 cosh 




| COSO) t - COSO) (t-t )| 

1 n n o. 1 


where the upper line corresponds to values t < t and the lower line, to the 

time t > t . 

o 


The free fluid surface displacement is 



av 




g 7T 

n=l 


( nir ^ 
co cos I — y 
n \ a J J 

r.2 


F(t) 


(3.3) 


where F(t) = sin co t + sin co (t - t Q ) ; the upper sign belongs to the time 

interval t < t while the lower sign belongs to t > t . 

o o 

The pressure in the container at a depth ( -x) is given by 
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At the tank wall y = 0 the pressure distribution yields 



(3.4) 


oo [~ n % “| 

a y n co n cosh L a (x + h)J 

7T 2 ’ 2 " 

n-1 rr cosh 


hr) 


F(t) 


and at the container wall x = a the pressure distribution is 


(3.5) 



-pgx - pV Q 



CO 



con cosh 


I HE. 

L a 


n 2 cosh 


r mr h Y 
i a / 


.(. x . t i 1 ). 


F(t) 


(3.6) 


At the container bottom x = -h the pressure distribution yields 
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^bottom = Pgh - Pv 0 


gh 



( WL \ 

\ a y; 


, a v / -i\ n cos 

+ tt 2 V 2 . /n7rh\ 

n=l rr cosh 

V a / 


F(t) 


(3.7) 


By integration of the appropriate pressure components, the liquid 
force and moment can be obtained. The force in y-direction becomes 


F = -mv 

y o 


2 v co _ , tanh 

— 2 2n-l 

n n=i (2n-l) 3 7 r h/ a 


(2n-d.br h 


F(t) 


(3.8) 


The fluid moment about an axis parallel to the z-axis through the 
center of gravity of the undisturbed liquid is 


mav °° co 


M z = 


7 r 


f z 


2n-i 


■L (2n - i) 


J 


n=i 


tanh 


|~ (2n-l)7r h j 


+ ? r 2 _ -j 

(2n-i)7rh/a [_ c °sh [(2n-l)7r h/ a] J 


F(t). 


(3.9) 


The velocity distribution in y-direction is given by 


8 $ 

u = — = v 
0y o 


I 

a 



+ if <-l) n cosh ^ sin 


n=l 


n cosh 


rary 


cosco (t -t) + cosco t 
no n 

cosco t - cosco (t-t ) 
n n o 


and in x-direction it is 
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u = 


3$ 

3x 


= V 


x+h 


o 


i V / \ n sinh (x + 

7r , (mrh\ 

n=l n cosh I— — I 


hj cos l™Lj\f C ° 301 n'V 0 + °° S “ / 


l a /< 


I COSO) t - COSOJ (t-t ) 
1 n n o 


B. RESPONSE OF LIQUID DUE TO A DOUBLE RECTANGULAR PULSE 


In considering a rectangular double pulse, i. e. , a pulse of con- 
stant magnitude v q for a duration of t seconds with a constant magnitude -v^ 

for a duration of t, seconds (Fig. 11) , 


v(t) = 


r 

V 

o 

for 

0 < t < t 

0 



p 

for 

t < t < t + ti 
o o 1 

) 


1° 

for all 

t > t + ti , 

o 1 j 


• 


This pulse can be described by the Fourier-Integral 


v(t) = 


(v + vA oo cos S3 1 sin £2 1 (v + v.\ oo sin S3 tfl -cos S2 1 \ 

^ 1 f 2_ ^ xU v_ r ' o) 


J 


a 


d£2 + 


/ 


Q. 


d£2 


oo cos £2 1 sin £2 ( t + L ) °° sin £2 1 Ti -cos 0 ( t + ti )1 

f 2 _l_d n-^f- L ° J 

7 T J £2 7 T J 

o o 


n 


d£2 


The velocity potential for this pulse is then: 
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■ llllllllll 


$ 




cos 12 1 sin 12 1 
c 

12 


d!2 


cos 12 1 sin 12 (t + ti) 
o 1 


12 


sin 12 1 £i-cos 12 tj 


_vi_ r fz__| _ -L (x 2 +2hx) ] 

7 r J L 2 a 6 2 a J 

o 

o 

oo _ . sin 12 1 f 1 - cos 12 ( t + ti)"| 

J- (x * + 2hx)] L- 2 — - J 

7 r ~ L2a 6 2a J 12 


d!2 


/ v + vi\ /, . 00 sin 12 1 cos 12 1 

l. P . ■ J rJ l ±\ 


IT 


S 


/ 


12 ' 


d!2 


( v o +V2 ) fh\ r°° Sinfit [ 1 _cosfit o] 

- — — e (-) J ^ ^ 


+ Mi) / 


00 cos 12 1 sin 12 ^t Q + 
3 


d!2 


12 ' 


d!2 


, x OC 

a / 


+ g 

7T 


^ 2a , , X 

+ p (v o + v,) 


sin Sit [l-cos Si (1 + tAJ 


12 ; 


d!2 


+ <V Vl) 

7 r u 


00 , fn 7 r , , , .1 (m. 'j °°12 cos 12t sin 12 1 

v> . ^ n cosh L"Bi(x + h)J cos \rt yj r o, 

A <_1 > t nir h J (co 2 -Si 2 ) df 

1=1 rr cosh o n 

a 

°° , Tux, , .1 (m. 'N 00 12 sin 12 1 fi-cos 12 1 I 

v' . c,n cosh L a (x + h)J cos V a y ) r L oJ 

h*- 1 * ; tut J (co 2 -si 2 ) 

n=l ir cosh — h o n 


2a y . ,.n cosh [a 2 " (x + h)l cos (’a'y) r ^ 

- 75 A*- 1 ) ; rnn ■> 

n=i rr cosh a o 


cos fit sin 12 (t + ti ) 
o 


(o> - 12 ) 

n 


dfi 


°° , Pn7r , , . 1 {nir \ °° fi sin fit fl-cos 12 (t + %)] 

. ?a V ,_!)■> cosh[T(x + h)JcosVr y ;|- l 2__ J dSi. 

7 T 3 1 2 , n 7 T, J (o> L ~ fi 2 ) 

n=i n cosh — h o n 

cl 
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In addition to the previous integrals others occur which can be solved 
by the same complex method indicated in the appendix. These are: 


/ 


°° cos fit sin fi (t + t[) 


fi 


dfi 


tt/ 2 for t < t + ti 
o 

for t > t + tt 


I 


/ 


sin fit £i-cos fi (t + W] 


fi 


'dfi 


sin fit £l-cos fi (t Q +ti)J 


fi ; 


dfi 


’ tt/ 2 for t < t + tj 


for t > t + ti 
o 1 


( 7r/ 4) t (2 t Q + 2tj - t) for t < t + tj 


(tt/ 4) (t +ti) ! 


for t > t + ti 
o 1 


/ 


°° fi cos fit sin fi (t + tj) 


a; 2 - fi 2 
n 


dfi 


(- 7 t/ 2 )cosoj tcosco (t +ti) for t < t + ti 
n no 1 o 1 

(tt / 2 ) sincu t sinco (t +ti) fort > t + t. 

n no 1 o ^ ■ 


/ 


sin fi (t + ti) cos fit 
o 

fi 3 


dfi 


(-tt/ 4) (^t 2 + (t* + t ) 2 ] for t < t + t t 


(~tt/ 2) t (t Q + t t ) 


for t > t + ti 
o 1 


r° fi sin fit |jL-cos fi(t D + tl} ] 

J t - » — * — 


w 2 - fi 2 
n 


dfi = 


(-7r/ 2) fcoso; t + sinco t sincu (t+ti)1 for t < t + b 

I L n n n o 1 J o 1 

( —TT / 2) ["COSCO t - COSCO t COSCO (t + till for t > t + ti 

^ L n n ,n o 1 J o 1 


With these results the velocity potential finally is found to be 
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“ 

#(y,x,t) = “ 

l 


= flL _ a _ J_ 


( x 2 + 2hx) v 


2a 6 2a ' J L ° 


Ml 

a 


+ vi 


— 


r 1 

i ( 2tt - 2t? - t 2 ) 
o 0 

V 

o 

< 

V 4 (t o - 2t) 



V 4 (t o - 2t > 

^ y 


r 

l 

4 

[(t o + t 1 ) 2 -2t 1 t-t Q 2 

< 

il 

i 

t 2 - 4 tt + 2t? - 2tfc 
“* o o J 



i[ t o 2+2ttl " (t o + tl)2 ] 


av , ,.n , fax, , (WL ) 

o y (-1) coshL a (x + h) J cos l a y) 

' * 2 nii n 2 oosh(S^) 

avi y ( -1) n cosh Da 1 (x + h)3 cos (a 2 ' y) 
^ n^l n 2 cosh ( yy " j 




COS w t + COS CO (t - t ) 

n n o 

cos co t - cos co (t - t ) \ 
n no/ 


cos co t - cos co (t - t ) 
n n o h 

v. J 




COS Cl) (t - t ) - COS Ct) ft - (t + ti )~ 1 

n o nL o 1 J 

-cos co (t - t ) - cos co ft - (t„ + t.j] \ 
n o n*~ ° liJ ( 

cos Ci) ft - (t + tj)”l - cos Ci) (t - t ) 
nL o J n o 


In this equation the first line corresponds to the time interval 0 < t < t Q , the 

second, to the time interval t Q < t < t + t* and the third line corresponds to 

t > t + ti . 
o 1 


In the symmetric case when the container wall at y = a is moved the 
same amount and with the same velocity to the left and right, the time inter- 
vals are consecutively 0 < t < t , t < t < 2t , and t > 2t . For this case 

J o o o o 

v q = vj, = t and the velocity potential is 
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t 2 - t 2 


t 2 - 4 t t + 3 t ^ 
o o 


o v ( 


-i) n coshCja?" (x + h)3 COS (¥v) 


n 2 cosh 


(3.10) 


cos co t + 2 cos co (t - t ) - cos co (t - 2 t ) 
n no no 


cos co t - 2 cos co (t - t ) - cos co (t - 2 t ) 
n no no 


cos co t - 2 cos co n (t - t Q ) + cos co n (t - 2 t Q ) 


The free fluid surface displacement is 


/ ^\ n+i 

v a °o (-1) co 


-2- < 2t -t > + 


(m l \ 
n C0S l a y j 


F (t) (3.11) 
n 


sin co t + 2 sin co (t - t ) - sin co (t - 2t ) 

n no n o i 

F (t) = < sin co t - 2 sin co (t - t ) - sin co (t - 2 t ) ) 

n i n no no 

sin co t - 2 sin co (t - t ) + sin co (t - 2 t ) 

V n no no-/ 


The upper line represents the time interval 0 < t < t^, the middle line, the 

interval t < t < 2 t and the lower line is valid for t°> 2t . Considering only 
o o o 
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the first term of the surface displacement, it can be seen that the surface level 
will sink during the time interval 0 < t < t and then rise again to its original 

level in the second time interval t < t < 2 t . The pressure in the container 

o o 

at a depth ( -x) is given by 



pv a ® 


(-D 


n 


co coshLa^(x + h 
n 


r lx , 


>3 


COS 


n 2 cosh 




F n (t) . 


(3.12) 


At the tank wall y = 0 the pressure distribution is 

r 

p = -Pgx-Pg(-) V ( 2 t - t 
y=o l a ] o ) o 


pv a °° (~l) n oj cosh Ij’a' ( x + h)] 

o y n 

7 r 2 9 , ( n7r h \ 

n=i n coshl— — I 


F (t) 
n 


(3.13) 


and at the container wall x = a it yields 

r 

t 


P y=a = - Pgx - pg 


a) % \ 


2 t - t 
o 


Pv a «> oj cosh Op <x + w ] 
o_ y n 

7 r 2 / - / , o , fair h 

n=l n cosh 


F (t) 
n 


(3.14) 


At tlie container bottom x = -h the pressure distribution is 


PbottonT PS h -MaJ V o 


\2 t - t 
o 


n / me ) 

pva « (-1) co cos ay 

_ o y n V / 

7T 2 A 2 , f n 7T h \ 

n=i n cosh | 

a ) 


V‘> 


The fluid force in y-direction is 
-2mv <*> co 


F = 
y 7T 


~ o y 2 n-l \. a 

2 h (2n-i) 3 7r (h/a) 

n=l 


tanh 


2n-i 


(t) 


( 3 . 15 ) 


( 3 . 16 ) 


and the moment of the liquid about an axis parallel to the z-axis through the 
center of gravity of the quiescent fluid is 


M z = 


mav «o co „ , 

o_ y 2n-l 

3 Lt (2n-l) 3 
7r 3 n=J v 


tanh 


^ (2n-l)7r h 


+ 


2 


(2n-l) 


7T h 

a 



( 3 . 17 ) 


The velocity distribution in y-direction yields 




— y 

7 r ^ 
n=l 


n 


m 


( -1) “ coshC^ (x + h)_ sin VC y> 

fmr h \ 
n cosh 1 

a ) 


G (t) 
n 


which vanishes for the left container wall y = 0; for the right container wall 

v - a it exhibits the value v during the time interval 0 < t < t and the value 
J o o 

-v during the time interval t < t < 2t , and vanishes for the time t > 2t , 
o 6 o o o 

as prescribed by the boundary conditions. 
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G n (t) represents the time function as given in the infinite series of 
the potential (3.10). The velocity distribution in x-direction is 



The first term in these expressions indicates the motion of the liquid without 
sloshing since the geometry of the container changes, The liquid has to flow 
down during the time interval 0 < t < t and flow up again during the time in- 
terval t Q < t < 2t Q as the wall moves back to its original position. 


C. RESPONSE OF LIQUID DUE TO SINUSOIDAL PULSE 

A sinusoidal pulse consisting of one single sine-wave (Fig. 11) 
of the form v^sin (27r t/ t Q ) , with t seconds duration as given by 


v(t) = jv o sin(2 7 r t/ t ) 


for 


o < t < t , 
o 

elsewhere 


can be represented by the Fourier-Integral 


2 v °° |"cos fit - 1 "1 cos fi 1 
o r L o J 


•<‘>=T "/ 


o o 


( S2 2 - 4 7T 2 / t Q 2 ) 


dS2 


2 v «> sin J2 t sin fit 
o 


-° / 


o o 

The potential is, therefore, 


(fi 2 - 4m 2 / t 2 ) 
o 


dfi . 
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— j-2- / flj - | - ^< x2 + 2hx > 

o 


J L 2 

O L 

COS fit p. -COS fit Q ^J 

(fi 2 - 47 r 2 / t Q 2 ) 


J 


sin fit sin fit 
o 

(fi 2 - 47T 2 / t 2 ) 
o 


2v gh oo sin fi t sin fit 


tifi - -r 1 — f 

t. a J 


t O a o' fi 2 (fi 2 - 47T 2 /t 2 ) 

o 


dfi 


2v gh oo cos fitT i-cos fit ~l 
o r L oJ 


t a 
o o 


/ 


fi 2 £fi 2 - 4 tt 2 / t Q 2 ] 


dfi 


^ aV o ^ . ,.n coshC^"(x + h)D Imt ) ? 

7T V' 1 ' 1 Jmh\ cos l aj 7 J 

o n=i n" 5 coshl — - — I o 


fi‘ 


cos fit£i-cos fit ~j 

(fi 2 - 4? r 2 / t 2 ) (w 2 - fi 2 ) 
o n 


dfi . 


fi 2 sin fit sin fit 
o 

(fi 2 - 4 t r 2 /t 2 ) (co 2 -fi 2 ) 
o n 


The occurring integrals are 


oo sin fit sin fit 


/ 


(0 2 - to Vt 2 ) an 


r 


o ■ 27T t 

T sin — 


"\ 


for t < t 


for t > t 


o ^ 


r 

cos fit f"l - cos fit "I , ~2T sin r 1 

°° L O- 1 ,^ l 4 t 0 

/ (fi 2 -47T 2 /t 2 ) 

^ o 


27 r t 


oo sin fit sin fit 


/ 


for t < t 


for t > t 


fi 2 / fi 2 - 4n 2 / t 2 
o | o 


dfi =< 


t Q I" • f ^7T t ^ - 7T t"| 

J 4rr 2 L 4 Sm ( to i 2 J 




for t < t 


o 

87T 


for t > t 
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fi 2 sin fit sin fit 


r 


(C0 n 2 - fi 2 ) / fi 2 - 47T 2 /t o 2 \ dn 

2 • ( 2^ A 

7T * Sin r~ t ^ 

-co n 77 r n 'o / 

0 , “2 — ^ — T7~L 2x Sin co t • cos co t + — 2 ^ — 27 — g v for t < t 

2(co - 4n V t l ) \ n no t (co -47 r Vt o 

n o L -J o v n o 


co 7 r 

2 n 2/ “2“ C0S w sin W Jr 


2(cu 2 - 47T 2 /t 2 ) 
n o 


n no 




for t > t 


cos fit [7 -cos fit 

fi 2 (n 2 _ 4 7 r 2 /t - 2 J > - dfi = 






f Z I 0 

, o . o . /2 tt if 

+ sr ( * ■ V - fir* sin rr 


for t < t 


for t > t 


o -J 


fi 2 cos fit[j.-cos fit J 

(cu 2 -fi 2 ) (fi 2 -% 2 /t 2 ) = 


r 


±7TC0 


n 


2(w 2 - 4 tt 2 /t 
n 


-r. sm co ' 
“) n 

o y L 


t + cos cu t sin cu t 
n no 


7i 2 sin t j 


for t < t 

2 ; 9 /. 9v O 


t (cu -47T 2 /t 2 ) 
o n o 


±7TCU 


n 


2 (cu 2 - 4 tt 2 /t 2 ) 
v. n o 


[ sin cu i 
n 


t T sin co t cos co t 
n no 


for t > t 



' With these results the velocity potential due to a sinusoidal pulse is 


4.(y,x,t)= v o [£- | - £ (*» + 2hx)] {j} sln^l 


sin 




2 v g t 2 h 
o o 


87T 


1 


2- + 1 


87 r 4nr t 


o 


87 r 


4a v °o 

x <-i) 

o n=l 


n coshC^ L (x + h)D cos y) ^ 
n 2 cosh "j 


G (t) 
n 


where 


G (t) = 
n 


(3.18) 


r 


_7rC ° n ... „ 2 7 r 2 smfe t ) _ . . ' 

2(co 2 -47r 2 /t 2 ) C Sm W n* + Sm " n ( t (co 2 - J/t *j f ° r 1 < *o 

i n o o n o 


7TC0 


2(co 


~ 2 — z — i > ~~2~ [-sm co t + sm co (t - t ) 1 
l - 477 V t *) L n n o -> 

n o 


for t > t 


The free fluid surface displacement is 



2v 7T 
° / 

x — - 

*‘o ' 


2v a 
0 


**'o 

where 

A 

e 

P Cs9 

—I 

F (t) = ( 

n 

L 

CO 2 

n 




y 2 _ a' 

2a 


XUL 

\a y 1 


cos 


'¥-‘H 

< O J J 


(3.19) 


n cos 


a 2 (co - 47 r 2 /t 2 ) * n 
n o 


F (t) 


cos 


27T t 


"N 


t + COS CO (t - t ) 

a n o 


] 


for t < t 


for t > t 


o 
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The pressure in the container at a depth (-x) is given by 



2pa(u.) 00 . n “I hzl 

__p y (-1) cosh L a (x + h) J cos \, a yj ^ 

7T t J 9 , /n L7T h \ , 9 , 9 / , n^‘ 

o n=i n J cosh (co -4irvty 

\ a. J n o 

At the tank wall y = 0 the pressure distribution is: 

2p V [~a I.?. . f ^ \ / 27T t \ 

— |j + 51 (3# + lo J cos (— j 


P y =o= -PS* + 


f r iij- _ iii 

pv ght J C0s t J 

o o \ <- o J 


2a7r 


2pv a oo 

~Tt~ 2 

o n=l 


n cosh ( x + h) 1 

9 , fm h ^ ( 9 47r 2 

n cosh^— j ^ n 2 -TT 


F (t) 
n 


and at the container wall y = a the pressure distribution yields 
2 PV 7T 


2pv 7r f “I f co 

pv ght 


cos ( 2tt t / t Q ) 


W - pw 


o o 
27r a 


2pv a °o 

- — — y 

7T t Li 


coshDa 1 (x + h)l 


o n=l if cosh 




2 / + 2 \ n 


F (t). 


(3.20) 


(3.21) 


(3.22) 
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The pressure at the bottom (x = -h) yields 


2 pv 7 r 


P bottom t 


7 r r- 

o 

< 

o L' 


_ a + h^_ 


2a 6 2a 

PV„ght, fZ COS ( 2 Tt/t 0 ) -1] 

0 


'"COS ( 27T t/ t Q ) ' 
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o y> (-1) cos 

7T t A 
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(3.23) 


The liquid force in y-direction is 
v v m 

F = --j ( 

y *o lo 


4v m °° I 42 n-al) 7 r h 

o ^ tanh V a / 


7T 


fc o n=i 7T f— ] (2n-i) 3 (co 2 -47r 2 /t 2 ) 
\,a/ 2 n-i o ' 


F 2„-l (t > 


(3.24) 


and the fluid moment yields 


M = ma 2 
z 


V 7T 
O 


COS (27T t/ t Q ) 


•ffl 


4y oo 


7 r at 


& 2 


F 2n-l (t) 


tanh 


f f(2n-l)7r h \ " A 

(2n-l)7r h ) v. cosh V a J J 
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+ (2n-i) 7r — 
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(3.25) 


Finally, the velocity distribution in y-direction is 
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It vanishes for the left wall y = 0 and exhibits for the right container wall y = a 
the expression v 


lu roil wan y 

° a loJ 


sin (2ir t/ t ) , as given in the boundary conditions. 


In x-direction it is 



4v 00 
+ T^ 2 


n 


(-1) sinh 


^ (x + h) 


n=i n cosh 


f mr h \ 

l a J 


cos 


rm 1 — . . 

— y) s n (t) 


which satisfies the boundary condition at the container bottom x = -h. 


SECTION IV. NUMERICAL EVALUATION AND CONCLUSIONS 

The behavior of the liquid caused by harmonic excitation of one con- 
tainer wall does not present any new effects compared to those of the excitation 
of the total container. At the natural frequencies co the free fluid surface 

displacement, fluid force and moment, etc., exhibit singularities. For this 
reason the number of graphs showing this type of excitation has been held to a 
minimum. However, it is interesting to note that in the case of harmonic ex- 
citation of one side wall the nodal line of the free fluid surface no longer re- 
mains at y = a/ 2, as in the case of excitation of the total container. For an 
excitation frequency Q, below first resonance, the nodal line is in the right half 
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of the free fluid surface, while above resonance it is in the left half of the free 
surface. Figure 12 exhibits the shape of the free fluid for excitation frequencies 
S2 = 0.9coj, 1.1 co i and 0. 9 uj. 


The character of the height of the free liquid surface is well presented 
by the formulas for double rectangular and sinusoidal pulses (Equations 3. 11 and 
3.19). The first one probably best describes the effect of the rocket firing. 

Both excitations, however, should yield a similar liquid response in spite of the 
different analytical expressions. The disturbance of one container wall acts 
like a line source emitting waves of all wavelengths and frequencies. The 
numerical results are given for a tank of length a = 250 feet and a liquid height 
h = 16 feet. The natural frequencies for this container are given in Table 1 for 
n = 1, 2 30. 

A rectangular double pulse with a duration of one second is considered. 

During the first half -second of the pulse the container wall moves to the right 

and the liquid moves down the moving wall, while during the second half of the 

pulse duration the container wall moves to the left and the liquid moves up, 

forming a valley in front of the wall. After the pulse is completed, the liquid 

is still moving toward the left, separating from the wall and starting to build 

up wave packages. As the first wave of the package moves towards the left, it 

decreases its height and the energy thus released is used to form other waves. 

The first wave valley of larger wave length initially created by the motion of 

the container wall moves faster away from the disturbance than the waves of 

smaller wave length in the wave package. The free fluid surface elevation of 

the liquid in the container is shown in Figures 13 through 24 at various times; 

the abscissa is (y/a) and the ordinate is x/ v . 

o 

It can be seen that after completion of the pulse the liquid reaches a 
height of x = 0. 41v q , while the valley exhibits a magnitude of x = -0. 042v q and 

is at a location y « 0. 95a. One second afterward, the valley has moved to 
y = 0. 87a and has decreased its magnitude to about x = -0. 038v q . At the same 

time, the peak amplitude of the wave which separated from the wall has moved 
only by 0. 02a and has decreased its amplitude to a value of x = 0. 212 v . As 

time increases, other waves are generated on the right side of the container, 
while the amplitudes of the present waves decrease continuously. 

At three seconds, i.e. , two seconds after the disturbing pulse has 
been completed, the original valley has a depth of x = 0. 028 v^ and is located 
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further to the left at y = 0. 77a. A second wave has been formed very close to 
the first one. The positive displacement of the first wave is reduced again to 
a value of approximately x = 0. 125 v , while its negative displacement is 

x = -0. 18 v . As time elapses, the waves move with different speeds toward 

the other container wall. The first valley needs about eleven seconds to traverse 
the container. 

At time four seconds, the first valley has further decreased its mag- 
nitude (to x = -0. 023v o > and is located at y = 0. 69a. The magnitude of the 

first wave has decreased to x - 0. 087 v q and is located at y = 0. 87a. A third 

wave has been formed and the second wave has a negative displacement of 
x = -0. 145a. As the first valley travels toward the left container wall its geo- 
metric shape becomes flatter and its depth decreases. The first waves separate 
more and more from the wave package, decreasing their amplitude, while others 
are still formed. 

At the time of six to seven seconds after completion of the pulse, the 
liquid on the left wall experiences a decrease in level. A few seconds later the 
first wave package seems to be completed. The first valley finally reaches the 
left wall at about eleven seconds with a magnitude of x =-0. 22 v , is reflected 

within two seconds and in the fifteenth second exhibits an amplitude of 
x = 0. 055 v . The reflected wave is running toward the right into the first waves 

of the wave package. 

At this time, about one -tenth of the container surface from the right 
hand wall is about at the equilibrium position x = 0. The first of the reflected 
waves increases in magnitude, then decreases during the time the waves work 
at each other. During that same time period, one-tenth of the liquid level on 
the right hand side container wall sinks slowly to a value of x = -0. 01 v^ at 

twenty-one seconds, and forms a curved surface with more negative displace- 
ment at the right container wall, while on the left container wall large liquid 
oscillations (x = -0. 06v q ) take place. With increasing time the liquid level on 

the right container wall rises again to a level of about x «= 0. 0l2v o , which 

stretches to about one-fifth of the fluid surface area from the right container 
wall into the liquid. 
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During the time period from twenty -five to twenty-seven seconds, 
the liquid surface rises at y = a but decreases at y = 0. 8a exhibiting an ampli- 
tude at y = a of about x = 0. 036 v , then swings back to the negative side to an 

amplitude of x = -0. 04 v , and so on. On the left hand wall of the container the 

amplitude increases to a value of about 0. 12 v q as time proceeds, while in the 

center of the container a "confused" wave motion takes place. 

At forty-nine seconds the amplitude of the liquid at the right container 
wall (y = a) exhibits a value x = 0. 068 v q and shows a wave of larger wave 

length in its vicinity. It should be noted, however, that for a fairly large value 
of v q = 10 feet per second, the original wave height at separation from the 

moving container wall would be 4. 1 feet, while the initial valley is only 0. 42 
foot in depth. The maximum wave height at the left container wall would be 
about 0.7 foot. 

Although the liquid motion can still be described for long periods 
after completion of the excitation pulse, in reality the motion may have damped 
out due to viscous forces of the fluid. Therefore, no further attempt is made 
to ascribe real significance to the wave motion which takes place after the first 
wave has reached the opposite container wall. 

Before proceeding to the description of the liquid behavior at a fixed 
location in the container, the results of the sinusoidal pulse excitation will be 
investigated and discussed. It should again be mentioned that no energy dissi- 
pation h as been considered during the liquid motion, indicating that the energy 
introduced into the system by the wall motion remains indefinitely in the liquid. 
This, of course, is not in agreement with daily experience. 

The response of the free fluid surface to a sinusoidal pulse of the 

form v sin ( 27 r t/ t ) of the container wall at y = a can also be presented at a 
o o 

fixed time after the pulse. Most of the results will be omitted here and only a 
few graphs are presented. (Figs. 25-30) . 

The total pulse duration is again t = i second. During half of the 

pulse duration the liquid will move down on the right container wall, while dur- 
ing the second half of the pulse, when the right container wall moves toward the 
left again to its initial position, the fluid moves up the wall. At time one and 
one-fourth seconds, i.e., one-fourth second after the pulse has been completed, 
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the valley exhibits a depth of x = -0.032 v q (which for v q = 10 feet per second 

is a valley of 0.32 foot) . The wave which moves away from the right hand wall 
has height of x =0. 25 v q (2. 5 feet for v q = 10 feet per second. 

A comparison of these results with those for the rectangular pulse 
shows that the valley and the first wave height exhibit about sixteen percent 
larger values in the rectangular pulse case. Essentially the same liquid be- 
havior is observed in these results. The first wave of the wave package moves 
with decreasing height but increasing wave length toward the left wall of the 
container, while new waves are added to the wave package as time elapses. 

The first created valley and waves move faster toward the left container wall 
than those of the interior of the wave package. 

At about six seconds the liquid level on the left container wall starts 
to sink until the first valley which decreased its original magnitude of 
x = -0. 032 v q to a value below -0. 01 v , increases it slightly again to about 

x = -0. 014 v at eleven seconds, at which time the reflection starts. Approx - 
o 

imately one-tenth of the liquid level on the right container wall is in the 
equilibrium position at that time. As the amplitude of the reflected wave on the 
left container wall increases to positive values of magnitude x = 0. 035 v , the 

liquid on the right hand side performs no motion at all from y = 0. 9a to y = a. 

As time progresses, the reflected waves work against the wave pack- 
age while the flat liquid surface area on the right hand side of the container 
increases its area from y = a to y « 0. 85a and sinks like a straight surface. 

In the center of the container a violent wave motion takes place. The wave height 
on the left container wall increases to x = ± 0. 04 v , while the right-hand side 

surface curves to larger depths for y = a, swings first up at y = 0. 8a and is 
then followed by an increase of surface elevation at y = a to a magnitude of about 
x = 0. 01 v q at twenty -five seconds. In the time proceeding, the free fluid sur- 
face elevation at y = 0. 8a decreases to negative values, while it increases at the 
right container wall to a magnitude of about 0. 024 v q at twenty-seven seconds. 

On the left container wall the amplitude is increasing with time and shows a 
value of x = 0. 07 v q at forty -two seconds. On the right container wall, a wave of 

larger wave length forms. 
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The sequence of these motions can be seen in a 16mm movie, which 
was produced from the computer results.' 1 ' It should be mentioned again, that 
the liquid motion is subjected to internal as well as to wall friction; therefore, 
the motion will be damped out soon and probably will not continue for the time 
length indicated in the graphs or the film. It should also be mentioned that the 
scale in the direction of the ordinate is greatly magnified. While the distance 
between the left and the right container wall is 250 feet, the height of the waves 
is only three to four feet. 

The flow rate of 4000 feet 3 per second during the firing of the rocket 
engines is simulated by a velocity v q of the right container wall of 


v = 1.67/ 1 . 
o o 

With this relation, the velocity amplitude v q of the right container wall during 
a double rectangular pulse for a half -pulse duration of t seconds has to be 


v = 6.7 
o 

ft/ sec. 

for 

t = 
o 

5 second 

v =3.4 
o 

ft/ sec. 

for 

t = 
0 

i second 

v = 1.67 
o 

ft/ sec. 

for 

II 

o 

+-> 

1 second. 


For a sinusoidal pulse, the value for v q can be obtained from 


16 ft x 150 ft. x 


ro ^ ^ f2n t\ 

J 


dt = 4000 ft 3 


(the width of the container is 150 feet) 


* 


The motion picture can be obtained on a loan basis from F. Kramer of 
MSFC, Test-Division. 



which yields 


1.67tt _ 5.25 

V o t t 

o o 

With this result the velocity amplitude v q of the container wall during a 

sinusoidal pulse of duration t seconds is obtained as 

o 


v = 10.5 ft/ sec for t = h second 
o o 


v = 5.25 ft/ sec for t = i second 

o o 


v = 2.63 ft/ sec for t = 2 seconds, 

o o 

For a rectangular -double pulse the original wave height at separation 
of the disturbing container wall would be x = 2. 74 feet, while the maximum 
liquid height at the left container wall is about 0. 48 foot for double-pulse dura- 
tion of one second. For a half-second duration, these values represent only 
the magnitude of 2.2 feet and 0.37 foot respectively. 

Similar results are obtained for a sinusoidal pulse; only those for the 
one-second duration are presented here. The liquid height at the right container 
wall is i.3 feet at the time of pulse completion, while the maximum wave height 
at the left container wall is about 0.37 foot. It may be mentioned that the 
sinusoidal excitation mode with one-second duration (t = 1 second) was used 

for the film. Figure 31 shows the surface elevation at the left wall, as a func- 
tion of time. 

The pressure distribution at the left container wall y = 0 is given for 
the rectangular pulse with a total duration of one-half, one and two seconds 
respectively, for various locations along that wall. The fluid height in this 
case is h 3 = 16 feet. For the pulse duration of one-half second, the pressure 
P + Pgz 

distribution -* 1 at the left wall y = 0 is given at the free fluid surface, 

pau G 
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at half the depth of the liquid and at the container bottom (Shown in Figs. 37 

through 39) . It can be seen that no dynamic pressure exists during the first 

few seconds that the wave is travelling from the right to the left container wall. 

At about nine seconds the pressure is indicated for the first time and continues 

to grow until it reaches (in a kind of sinusoidal fashion) a maximum of about 

0. 0035 p a v at the free fluid surface, 0. 0014 p a v + p g 0. 032a at half the 
o o 

liquid depth and 0. 0013 p v q a + p g 0. 064a at the container bottom. 

However, these peak pressures are reached at various times, as 
seen in Figures 32-40 for various pulse durations. With the mass density p * 0.973 
lbs sec 2 / ft 4 and the specific weight of y = p g = 62. 4 lbs/ ft 3 , the total peak 
pressures at the left container wall for a rectangular double pulse are 


(a) for 

V * 

second by 


11. 

7 lbs/ 

ft 2 

at 

X = 

0 

504. 

7 lbs/ 

ft 2 

at 

X = 

-ih 3 

1004. 

7 lbs/ 

ft 2 

at 

X = 

-h 3 

(b) for 

V 3 

r second by 


21 

lbs/ 

ft 2 

at 

X = 

0 

511. 

7 lbs/ 

ft 2 

at 

X = 

" h 3/ 2 

1008. 

4 lbs/ 

ft 2 

at 

X = 

-h 3 

( c) for 

t =1 
0 

second 

by 


27 

lbs/ 

ft 2 

at 

x = 

0 

518. 

5 lbs/ 

ft 2 

at 

X = 

" h 3/ 2 

1017 

lbs/ 

ft 2 

at 

X = 

-h 3 


(at the free fluid surface) 
(at half the depth) 

(at the container bottom) 


(at the free fluid surface) 

( at half the depth) 

(at the container bottom) 


(at the free fluid surface) 
( at half the depth) 

(at the container bottom) 


It can be seen, that the pressure at the bottom of the container (Fig. 
41) coincides with the pressure at the wall (y = 0 for x = -h) ( Fig. 39) . In the 


center of the container at y * a/2, the pressure ( 


^Bottom 

av p 
o 


) versus time 


is presented in Figures 44, 45, and 48 for t = 1/ 4, 1/ 2 and 1 second. The 
peak pressure at the container bottom at y = a/ 2 is, therefore, 
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1002.7 lbs/ ft 2 for t = 1/ 4 second 

o 

1005.2 lbs/ ft 2 for t = 1 / 2 second 

o 

1009.6 lbs/ ft? for t = 1 second . 

o 

The fluid force F^ is exhibited in the Figures 50 through 52 for 

t = \ and 1 second. For t Q = J second, the force per foot of tank width in 

y-direction drops from a value of 113 pounds per foot to 61 pounds per foot. 

For a tank width of 150 feet this yields a total force on the container in y-direc- 
tion of 17000 pounds at t = \ sec. , i. e. , exactly at the end of the pulse, and 
9100 pounds at about 48. 5 seconds. For a half -pulse duration t = \ second, 

the force per foot drops from 205 pounds per feet at t = 1 second to 115 pounds 
per feet at about 48.5 seconds. This corresponds to a total force 30, 800 pounds 
and 17,200 pounds, respectively. For t = 1 second, i.e., a pulse duration of 

two seconds, the force per foot yields 263 pounds per foot at t = 2 seconds and 
drops to 182 pounds per foot at about 49 seconds. This corresponds to a total 
force of 39, 400 pounds and 27, 300 pounds, respectively. 

The liquid moment is presented in Figures 53 through 55. For 
t = 5 second, the maximum value of the moment per foot of liquid is 6700 

pounds, which corresponds to a total moment of 10 6 foot-pounds. For t = \ 

second the maximum value is 14, 300 pounds corresponding to a total moment 

about the center of gravity of the undisturbed liquid of 2, 142, 000 foot-pounds. 

For t =1 second the maximum value of the moment per unit width of the con- 
o 

tainer is 26, 300 pounds. This corresponds to a total moment of M =3, 950, 000 

z 

foot-pounds . 

SECTION V. RECOMMENDATIONS 

According to the theory presented and its numerical evaluation, it 
can be recommended that the walls of the container should be at least three 
feet higher than the liquid height to prevent any liquid loss during the operation 
of the system. The wall structure of the liquid reservoir should be designed 
according to the maximum pressures listed on page 35. 
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Since it was found that the liquid force and moments for pulses 
caused by one -wall -excitation are considerably smaller than those of a har- 
monic excitation of the complete container, further design specifications can 
be found only by treating the total system of liquid and structure. This has 
been done for various wind inputs and the result is presented in the report 
"Interaction of Structure and Liquid in the Sound Suppressor System. ” 


George C. Marshall Space Flight Center, 

National Aeronautics and Space Administration 
Huntsville, Alabama, September 30, 1965. 


37 



TABLE I. NATURAL CIRCULAR FREQUENCIES 




y-Direction 



z -Direction 

0 0 

n 

hi / ai 

h>/ a.i 

ha/ a! 

hi/ a. 

h 2 / a 2 

ha/ a 2 

fluid height ratio 

COi 

0.2248 

0.2558 

0.2832 

0.4293 

0.4863 

0.5319 

co 2 

0,4461 

0.5051 

0.5555 

0. 8354 

0. 9305 

1. 0044 


CJ 3 

0.6608 

0. 7422 

0.8087 

1.2024 

1.3104 

1.3837 


0>4 

0. 8663 

0. 9634 

1.0380 

1 . 5243 

1.6358 

1.6851 


“5 

1.0609 

1.1666 

1.2425 

1 . 8028 

1. 8880 

1. 9306 



1.2435 

1.3515 

1.4236 

2.0439 

2. 1104 

t 

2.1389 


CO 7 

1.4136 

1.5191 

1. 5842 

2.2548 

2.3042 

2.3223 


co$ 

1.5716 

1.6710 

1.7276 

2 . 4421 

2.4774 

2.4886 


CO 9 

1.7180 

1. 8090 

1. 8569 

2.6111 

2.6358 

2.6425 


10 

1. 8538 

1. 9353 

1 . 9748 

2.7660 

2.7829 

2.7869 


Wi 1 

1. 9797 

2. 0515 

2.0835 

2. 9098 

2. 9213 

2. 9236 


w 12 

2.0970 

2.1592 

2.1848 

3.0450 

3.0526 

3.0539 


WiS 

2.2066 

2.2598 

2.2801 

3.1730 

3.1780 

3.1788 


«14 

2.3094 

2.3545 

2.3704 

3.2951 

3.2984 

3.2988 


w 15 

2.4062 

2.4441 

2.4565 

_ 3.4123 

3.4144 

3.4147 


Wie 

2.4979 

2.5295 

2.5390 

3.5251 

3.5265 

3.5267 


CO 17 

2.5851 

2.6113 

2.6186 

3.6342 

3.6351 

3.6352 


w 18 

2.6683 

2.6899 

2.6955 

3 . 7400 

3.7406 

3.7406 ’ 


w 19 

2.7480 

2.7657 

2.7700 

3.8427 

3. 8431 

3.8431 


w 20 

2.8247 

2.8392 

2. 8424 

3.9427 

3. 9429 

3. 9429 


W 21 

2. 8987 

2. 9105 

2. 9129 

4. 0402 

4. 0403 

4. 0403 


w 22 

2.9073 

2. 9798 

2. 9817 

4. 1353 

4.1354 

4.1354 


w 23 

3.0397 

3.0474 

3.0488 

4. 2283 

4.2283 

4.2283 


^24 

3. 1072 

3.1135 

3.1145 

4.3192 

4.3193 

4.3193 
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TABLE I. (Concluded) 




y-Direction 



z -Direction 

05 

n 

h l/ a i . . 

h _2/ a i 

h- 3 /a. 1 

h l/ a 2 _ 

h 2 /a 2 

h 3 /a 2 fluid height ratio 

w 25 

3.1730 

3.1780 

3.1788 

4. 4083 

4. 4084 

4. 4084 

w 26 

3.2372 

3.2412 

3.2418 

4.4956 

4. 4957 

4. 4957 

05 2 7 

3.2999 

3.3031 

3.3036 

4. 5813 

4.5813 

4.5813 

^28 

3.3613 

3.3639 

3.3642 

4. 6654 

4. 6654 

4. 6654 

0)29 

3.4214 

3.4235 

3.4238 

4. 7479 

4. 7479 

4. 7479 

w 30 

3.4805 

3.4821 

3.4823 

4. 8291 

4.8291 

4. 8291 



h* = 10 

feet 

a l 

= 250 feet 



h 2 = 12 

feet 

a 2 

= 150 feet 



h 3 = 16 

feet 
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FIGURE 2. TANK GEOMETRY AND COORDINATE SYSTEM 
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FIGURE 9. PRESSURE 
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FIGURE 11. PULSE EXCITATIONS 
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FIGURE 13. LIQUID SURFACE ELEVATION FOR RECTANGULAR DOUBLE 
PULSE OF ONE SECOND DURATION FOR VARIOUS TIMES 
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FIGURE 14. LIQUID SURFACE ELEVATION FOR RECTANGULAR DOUBLE 
PULSE OF ONE SECOND DURATION FOR VARIOUS TIMES 
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FIGURE 15. LIQUID SURFACE ELEVATION FOR RECTANGULAR DOUBLE 
PULSE OF ONE SECOND DURATION FOR VARIOUS TIMES 
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FIGURE 16. LIQUID SURFACE ELEVATION FOR RECTANGULAR DOUBLE 
PULSE OF ONE SECOND DURATION FOR VARIOUS TIMES 
















FIGURE 18. LIQUID SURFACE ELEVATION FOR RECTANGULAR DOUBLE 
PULSE OF ONE SECOND DURATION FOR VARIOUS TIMES 
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FIGURE 19. LIQUID SURFACE ELEVATION FOR RECTANGULAR DOUBLE 
PULSE OF ONE SECOND DURATION FOR VARIOUS TIMES 
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FIGURE 20. LIQUID SURFACE ELEVATION FOR RECTANGULAR DOUBLE 
PULSE OF ONE SECOND DURATION FOR VARIOUS TIMES 


59 









6 







o 



y/ 


FIGURE 22. LIQUID SURFACE ELEVATION FOR RECTANGULAR DOUBLE 
PULSE OF ONE SECOND DURATION FOR VARIOUS TIMES 
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FIGURE 23. LIQUID SURFACE ELEVATION FOR RECTANGULAR DOUBLE 
PULSE OF ONE SECOND DURATION FOR VARIOUS TIMES 
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FIGURE 24. LIQUID SURFACE ELEVATION FOR RECTANGULAR DOUBLE 
PULSE OF ONE SECOND DURATION FOR VARIOUS TIMES 
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FIGURE 25. LIQUID SURFACE ELEVATION FOR SINUSOIDAL PULSE OF 
ONE SECOND DURATION FOR VARIOUS TIMES 
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FIGURE 26. LIQUID SURFACE ELEVATION FOR SINUSOIDAL PULSE OF 
ONE SECOND DURATION FOR VARIOUS TIMES 
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FIGURE 27. LIQUID SURFACE ELEVATION FOR SINUSOIDAL PULSE OF 
ONE SECOND DURATION FOR VARIOUS TIMES 
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FIGURE 28. LIQUID SURFACE ELEVATION FOR SINUSOIDAL PULSE OF 
ONE SECOND DURATION FOR VARIOUS TIMES 
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FIGURE 29. LIQUID SURFACE ELEVATION FOR SINUSOIDAL PULSE OF 
ONE SECOND DURATION FOR VARIOUS TIMES 
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FIGURE 31. LIQUID ELEVATION VERSUS TIME AT y = 0 FOR RECTANGULAR 
DOUBLE PULSE 
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FIGURE 33. PRESSURE DISTRIBUTION AT THE LEFT CONTAINER WALL FOR 
RECTANGULAR DOUBLE PULSE FOR VARIOUS PULSE DURATIONS 
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FIGURE 34. PRESSURE DISTRIBUTION AT THE LEFT CONTAINER WALL FOR 
RECTANGULAR DOUBLE PULSE FOR VARIOUS PULSE DURATIONS 
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FOR VARIOUS PULSE DURATIONS 
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FIGURE 36. PRESSURE DISTRIBUTION AT THE LEF’ 
RECTANGULAR DOUBLE PULSE FOR V, 
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FIGURE 37. PRESSURE DISTRIBUTION AT THE LEFT CONTAINER WALL FOR 
RECTANGULAR DOUBLE PULSE FOR VARIOUS PULSE DURATIONS 
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FIGURE 41. PRESSURE DISTRIBUTION AT THE CONTAINER BOTTOM FOR 
RECTANGULAR DOUBLE PULSE 
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FIGURE 42. PRESSURE DISTRIBUTION AT THE CONTAINER BOTTOM FOR 
RECTANGULAR DOUBLE PULSE 
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FIGURE 44. PRESSURE DISTRIBUTION AT THE CONTAINER BOTTOM FOR 
RECTANGULAR DOUBLE PULSE 
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FIGURE 46. PRESSURE DISTRIBUTION AT THE CONTAINER BOTTOM FOR 
RECTANGULAR DOUBLE PULSE 
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FIGURE 47. PRESSURE DISTRIBUTION AT THE CONTAINER BOTTOM FOR 
RECTANGULAR DOUBLE PULSE 
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FIGURE 48. PRESSURE DISTRIBUTION AT THE CONTAINER BOTTOM FOR 
RECTANGULAR DOUBLE PULSE 
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FIGURE 50. LIQUID FORCE FOR RECTAN 
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FIGURE 51. LIQUID FORCE FOR RECTANGULAR DOUBLE PULSE FOR 
VARIOUS PULSE DURATIONS 
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FIGURE 52. LIQUID FORCE FOR RECTANGULAR DOUBLE PULSE FOR 
VARIOUS PULSE DURATIONS 
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FIGURE 53. LIQUID MOMENT FOR RECTANGULAR DOUBLE PULSE FOR 
VARIOUS PULSE DURATIONS 


92 




FIGURE 54. LIQUID MOMENT FOR RECTANGULAR DOUBLE PULSE FOR 
VARIOUS PULSE DURATIONS 
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FIGURE 55. LIQUID MOMENT FOR RECTANGULAR DOUBLE PULSE FOR 
VARIOUS PULSE DURATIONS 


94 




vO 

Ln 


FIGURE 56. PATH OF INTEGRATION FOR IMPROPER INTEGRALS 




APPENDIX 


EVALUATION OF IMPROPER INTEGRALS 

The integrals in the text are improper integrals and require the 
evaluation of 


/f(S2) e iaS2 dfi 
C with a > 0. 

The path C of integration is the real axis and the semicircle around the origin 
in the upper half plane ( Fig. 56) . 

Evaluation is restricted to two of these integrals. The others are very similar 
and do not require additional techniques. 

The integral 


oo 


/ 


cos £2t sin £2 1 
£2 


d£2 


2 


I 


O 


sin .g - (^T J .^ 
a 




o 


because sin a T cos Q t = 1/ 2 sin £2 (t + T) + 1/ 2 sin £2 (T - t) . Essentially, 
it is necessary to determine the value of the integral 


i 


o 


sin a £2 

a 


d£2 


since a can be taken to be either T+t or T-t. 


I 


sin a Q 


dsJ - i I 


+°° iafi 




dS2 - 


It is 

+oo — i a Q. 



oo 
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The integral 


I 


£2 cos £2t sin £2 T 
u 2 - £2 


i r £2 sin £2 (t+T) 

dQ " * J co 2 - n 2 

o 


d£2 + 


yields essentially integrals of the form 

H-oo 


/ s 


-oo co - £2 


ia£2 

2 •' 


1 r o sin Q ( T -t) HO 

2 J co 2 - a 2- 


These integrals exhibit singularities on the integration path but their 
principal value can be found by complex evaluation methods. The integration 
path is taken along the real axis and a semicircle in the complex £2 -plane. If 
it can be shown that the integral along the semicircle vanishes with increasing 
radius of the circle, the value of the integral along the closed path is the same 
as that of the integral along the real axis from - » to +°° . If the integral 


f f (£2) e lfia dfi 
Ci 

is considered with a > 0 and Cj a semicircle around the origin in the upper half 
plane, then f ( £2 ) may satisfy the following conditions : 

1. f( £2 ) be analytic for Im( £2 ) > 0 except for a finite number of poles. 

2. The absolute value of the function f (£2) on the curve approaches zero 
as the radius of the circle C* approaches infinity. 

lim j f( £2 ) | = 0 for 0 < arg £2 ^ F. 

| 1 ^ <=o 

Then 

lim j'f(tt) e ifia d£2 = 0 , if a > 0 . 

|S7 [ — ► oo Ci 
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Proof : With the second condition, it can be concluded with sufficiently large 
|£2| , that |f(£2)| < e for all £2 on Ci and e > 0. 

With £2 = Re * ^ it is 

I i£2a| I ia R(cos cp + i sin cp) I _ -aR sin <p -2a Rcp/ir 

\ B I |0 | — 0 ^ G 

because of 


2 / 7 r s < 1 for 0 < cp < 7r/ 2. 

<P 

Therefore, it is 


f f( £2 ) e lS2a d£2 

= 

7T 

f (/n 1 a R cosw -4 R sin cp „„ i<p , 

J f ( Re ) e .e ^.R£2e^d<p 

,Ci 


0 


7T „ . 7r/2 _ . 

. -a r -a R sin w r - a R since* , 

< eR j e r dg> = 2 eR j e ^ dcp 

o 0 


_ „ -of -2a J 7T£ . -a R 

< 2<R J e d^> = — 1-e 


The limit value for R — °o ( e — 0) yields 


lim f f(£2) e iaS2 d£2 

R — - 00 Ci 


0. 


For the type of integrals appearing in the text the function f(£2) is a 
rational function, consisting of two polynomials h( £2 ) andg(£2). It is 


f(Q) 


= h(£2 ) 
g(S2) 


99 


with at least the order of g( £2 ) being larger by one than that of h( £2 ) , i. e. , 



Therefore, 

lim |f(£2)| = 0. 



3 , 

If the function f ( £2 ) e satisfies the following conditions : 
ia £2 

1 . f ( £2 ) e is analytic in the upper half plane for Lm£2 > 0 except for 
a finite number of poles. 

Iq a 

2. f(£2) e has only simple poles on the real axis. 

3. lim | £2 f(£2) e* a ^ | = 0 for 0 £ arg £2 £ n 
£2 -*■ °° 

then 


f f(£2) e ia d£2 = 27r i ^ Res. (a ) + in ^ Res. (b ) 

A=1 X=1 

where the a 's are n poles in the upper half plane Im£2 > 0 and the b 's are m 
A . _ A 

iS2 cl 

simple poles on the real axis for the function f ( £2 ) e 

Proof: The first part of this statement follows from Cauchy's Theorem and the 
fact that the integral along the semicircle vanishes. Let us proof the second 
part of the formula. At £2 = ± co are simple poles on the real axis. A semi- 
circle of radius p is taken around £2 = ± w in the upper half plane. The contour 
C consists of the real axis from -R to +R ( excluding the segments 
±co -p to ±co +p) , the semicircles p and the semicircle in the upper half 
plane having the radius R and the center at the origin. With Cauchy's Theorem, 
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$f(S2) e iaS2 d£2 = *ff(C2) e dfi + /fit f (fi ) dfi + I 

C C t (Pi) n 


+ f(S2) e lS2a dn 

(Pi) 

n 

= 2 7r i X Res (a ) 

X=i 71 

where a are the poles inside the integration path C. 
A 

The integral I is 


I r = / 1 f(fi) e mS2 dS2 + / P2 f(C2) e ian dtt + f f(S2) e ia$2 

-R -w -Pi co+Pz 

Where P denotes the prinicpal value, 

+0O 

lim I = P f f ( £2 ) e ia ° dS2 
n J 

— oo 

R — ► 

Pi — 0 

In considering the integral 


■j\ f(£2) e* 3 ^ dfi 

(Pi) 

where, on the semicircle, can be written as 

+ £2 = . (jo + Pz ( 0 < cp < 7r ). 


d£2 . 
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let 

Since f (fi ) e has at fi = w a simple pole, the Laurent expansion of the func- 
tion can be written as 

oo 

„ . _ v iafi bj , v x M 

f(fi) e = 4 + 2 j 0 ! (i2 -co ) . 

w p = 0 M 

The residue at fi = c o is b 2 . Therefore, 

£f«l) e iaS2 d!2 = /\ i 4, + f / « ip? +1 e i(fl + 

(P 2 ) 7T M = 0 7T 

which yields with p 2 -*■ 0 

lim f(fi) e ia ^ dfi} = — i 7T b 2 . 

P2 -* 0 (P2) 

The same is valid for other simple poles on the real axis and proves the above 
statement. 


The preceeding integrals satisfy the cited conditions, exhibit only 
simple singularities on the real axis and are 


+00 

I 


— OO 


iafi 
e 

fi 


dfi 


7T i 


(a >0) 


and 


+00 

I 


— 00 


iafi 

fie 

3 • y " 
co fi 


dfi 



-iaw 


co e 

2co 


) 


— it i cos aw; (a > 0) . 


* 

For a < 0 the integration path C has to be taken in the lower half plane 
( 7 r < arg fi < 2 7r ) . Integration in clockwise direction then yields a negative 
sign in Cauchy's expression. 


102 



Now the first integral 


/ 


cos ft t sin £2 T 
dft 


dn = iT / 


, +<*> i(t+T ) ft . +■*> -i(t+T) n 

ire i r e 

■ dft - — J — dft 


ft 


ft 


+oo i(T-t) ft , +°° -i(T-t) ft 

f 2 dn .1 f ' 

8i J .O. 8i J 


ft 


dft 


which yields for t < T 


OO 


/ 


cos ft t sin ft T 
ft 


dft 


£ 

2 


For t = T the third and fourth integral vanish together and the value of the inte- 
gral is 7 r/ 4. For t > T the integrals cancel each other. It is, therefore, 


/ 


cos ftt sin ftt 
ft 


dft = 


V 

for 

t 

< 

T 

7T 

4 

for 

t 


T 

^0 

for 

t 

> 

T 


The second integral to be evaluated yields 


/ 


ft cos ft t sin ft T 
co 2 - ft* 


= - f 

8i J 


+°° i(t+T) ft 


fte 


CO 2 - ft 2 


dft - 


- ~ f 

8i J 


+°° -i(t+T) ft 


fte 


+=o i(T-t)ft 


CO 2 - ft 2 


. i r fte‘' A „ 

dft + TT I 2 ^T5 — ft 

8i J co - ft* 1 


8i 


+oo -i (T-t) ft 

r fte v ln 

9 ^2 dft . 

J 00 1 - ft^ 
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This is for t < T 


r £2 cos £2 tsin£ 2T w 7r ir 

I § 772 d£2 = cos co(t+T) - — cos co (t+T) - -cosco(t-T) 

Q CO “ o O O 

- ^ cos co (t-T) - - ~ cos co t cos co T. 

o Z 

For t = T the third and fourth integral vanish together and the value of the inte- 
gral is - it/ 4 cos 2 co t, while for t > T the exponential function exhibits in the 
third and fourth integral a negative value for (T-t) in the exponent, which results 
in a change of the integral path from the upper to the lower half plane and vise 
versa for the fourth integral. It is, therefore. 


/ 


£2 cos £2 t sin n T 7r n /jL im . , n ,, 

d£2 = - — cos co (t+T) - — cos co (t+T) + — cosco (t-T) 

o o o 


co 2 - £2 2 


+ 7 cos co ( t-T) = f - sin co t sin co T for t > T 

O A 


It is, therefore. 


I 


cos £2 t sin £2 T 
£2 


r. 


d£2 = 


<7 


for t < T 

for t = T ) 

for t > T 


/ 


£2 cos £2 t sin £2 T 
co 2 - £2 2 


d£2 = 


f 

£ 




- 

COS CO t COS CO T 

for 

t < T 


2 




<- 

7T 

4 

cos 2 cot 

for 

t =T 


£ 

2 

sin co t sinco T 

for 

t > T 


The other integral expression can be evaluated in a very similar way (see also 
Ref. 6). 
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edge of phenomena in the atmosphere and space. The Administration 
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